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ANGULAR  RESOLUTION  OF  COHERENT  AND 
NONCOHERENT  SOURCES 


INTRODUCTION 

When  attempting  to  resolve  targets,  one  can  either  resolve  them  in  angle  or  range  or  a  combina¬ 
tion  of  both  coordinates.  However,  the  range  coordinate  is  a  much  more  powerful  discriminant  for 
separating  targets.  For  example,  two  targets  separated  by  1000  ft  are  separated  by  two  l-/is  pulsewidths 
but  by  only  1/20  of  a  2°  beamwidth  at  a  range  of  100  nmi.  On  the  other  hand,  when  attempting  to 
resolve  radiating  sources  (e.g.  jammers),  one  must  resolve  them  in  angle. 

The  probability  of  resolving  sources  is  not  only  a  function  of  their  separation  but  also  a  function 
of  their  strengths  and  phase  differences.  There  have  been  many  articles  written  on  resolution  and 
reprints  of  some  of  the  most  important  ones  can  be  found  in  Ref.  1 .  In  an  early  article  Ksienski  and 
McGhee  [2],  using  a  decision  theoretic  approach,  indicate  that  targets  within  a  quarter  of  a  beamwidth 
can  be  resolved.  Lately,  the  problem  of  angular  resolution  has  been  investigated  by  using  superresolu¬ 
tion  techniques  [3].  All  of  these  investigations  solve  the  problem  by  using  estimation  approaches.  We 
approach  the  resolution  problem  as  a  binary-hypothesis  problem  where  the  two  hypothesis  are: 

•  H\.  One  Source  Present 

•  H2:  Two  Sources  Present 

The  development  is  similar  to  the  one  used  for  the  range  resolution  of  targets  [4]. 

GENERALIZED  LIKELIHOOD  RATIO  TEST — NONCOHERENT  SOURCES 

We  now  formulate  the  resolution  problem  as  a  binary-hypothesis  test.  The  received  samples, 
from  the  two  sources  impinging  on  the  array  elements  shown  in  Fig.  1,  are 

xu  “  %  +  -4/1  cos  ($/j  +  ir(i  -  1)  sin  flj) 

+  An  cos  (<£y2  +  vU  -  1)  sin  92), 

and  (1) 

Yu  —  nru  +  Aj\  sin  +  rrO  -  1)  sin  9]) 

+  Aji  sin  (<f> ji  +  ir(i  -  1)  sin  92 ), 

where 

•  Xy  and  YtJ  are  the  inphase  and  quadrature  samples  from  the  /th  element  and  yth  range 
cell, 

•  nX(j  and  nYfj  are  independent  inphase  and  quadrature  Gaussian  noise  samples  with  mean 
zero  and  variance  o-2, 

•  Aj\  and  Aj2  are  source  amplitudes  in  the  jih  range  cell, 

•  <t>j\  and  4>j2  are  the  phase  angles  in  the  yth  range  cell, 

•  9\  and  92  are  the  angles  of  arrival,  and 

•  the  receiver  elements  are  1/2  wavelength  apart. 


Manuacript  approved  March  7,  1984. 
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SOURCE  #1 


SOURCE  #2 


Fig.  1  —  Signals  from  two  sources  received 
by  an  8-eiement  linear  array 


O  O  O  0*0  o  o  o 


ARRAY  ELEMENTS 


The  angular  resolution  problem  is  equivalent  to  deciding  which  hypothesis  is  true: 

H\.  Aj\  >  0  and  Aj2  —  0 
or 

H2:  Aj\  >  0  and  AJ2  >  0. 

If  [Aji,  Aj2,  4* j\>  4>ji.  9\.  02)  are  all  known,  the  optimal  detector  (in  the  Neyman  Pearson  sense)  is  the 
likelihood  ratio.  Since  in  the  problem  of  interest  the  parameters  { AJt ,  Aj2,  <t>j2,  9it  02)  are  un¬ 
known,  in  the  likelihood  ratio  we  replace  the  parameters  with  their  maximum  likelihood  estimates. 

The  maximum-likelihood  estimates  are  the  values  of  the  parameters  which  minimize  the  square 

error 

m  n 

L  —  ( X(j  -  Aji  cos  (0yi  +  tt(/  —  1)  sin  9 j)  —  AJ2  cos  +  ir  (i  -  1)  sin  92)r 

+ls  (Y(j  -  Aji  sin  (^71  +  ttO  -  1)  sin  0i)  -  AJ2  sin  (<f>j2  +  ir(i  -  1)  sin  02))2,  (2) 

where  n  is  the  number  of  array  elements  and  m  is  the  number  of  range  cells.  Making  the  following 

substitutions; 

Bj  ™  Aji  COS  0yl< 

Cj  "  Aj2  sin  <f>ji, 

Wt  -  cos  (w(i  -  1)  sin  9\), 

Ui  «  sin  (ir(/  -  1)  sin  9 1), 

Dj  —  Aj2  cos  9j 2, 

Ej-  A j2  sin  9j2, 

St  —  cos  (ir(<  -  1)  sin  92), 


Tt  —  sin  (w (/  -  1)  sin  92), 
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Eq.  (2)  can  be  written  as 

L  -  2j  ij  (*</  -  BjW,  +  CjU,  -  Pj  S,  +  Ej  T,)2 

+  Jj  Jj  OV  -  C;  ^  -  Bj U,  -  £,S,  -  D}  T,)2.  (3) 

For  any  given  value  of  0i  and  IF,,  £/„  Sh  and  T,  are  known,  and  the  maximum-likelihood  esti¬ 
mates  of  Bj,  Cj,  Dj,  and  Ej  can  be  found  by  solving  the  equations 

BL  B  L  BL  BL  n  . 

BBj  "  BCj  “  BDj~BEj~  J  m- 

The  Appendix  of  this  report  shows  that  the  solutions  for  the  one  source  case  (that  is,  assuming  f/j  is 
true)  are 

Dj~  Ej-  0, 


Bj  -  Jj  OF,*„  +  Ut  Yu)tn, 


Cj-j jOF,  »</-  ATyVn. 


In  the  Appendix  it  is  also  shown  that  the  solutions  for  the  two-source  case  (that  is,  assuming  H2  is 
true)  are 

Bj  -  ( nHj  ~  FJj  ~  GKj)/(n2  -  F*  -  G2). 

Cj  -  (n/j  +  GJj  -  FKj)/(n 2  -  F2  -  G2), 

Dj  -  (nJj  -  FHj  +  GIj)/(n2  -  F2  -  G2),  (5) 


Ej  -  ( nKj  -  GHj  -  FIj)/(n 2  -  F2-  G2). 

where 

F-j^iW'S'  +  U'T,). 

G  -  Jj  (U,S,  -  W,T,). 
ffj  -  Jj  (W,XV  +  U,YU). 

/,-  U,XU), 

Jj-j^iS'Xu+T'Y"), 

and 

Kj  -  ij  iS,Yu  -  W- 

A  direct  search  technique  (51  is  used  to  estimate  the  source  directions  0|  and  92. 
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Finally,  the  log  likelihood  ratio  can  be  written  as 

A-  Min  Ul-  Min  {£}.  (7) 

Dj-Ej-0 

Calculation  of  the  log  likelihood  involves  two  minimizations.  For  the  first,  one  searches  on  6]  and  uses 
Eq.  (4)  for  the  other  parameters.  For  the  second,  one  searches  on  0\  and  02  and  uses  Eq.  (S)  for  the 
other  parameters.  Previous  work  (4]  on  range  resolution  indicate  that  results  similar  to  those  obtained 
with  Eq.  (7)  can  be  obtained  by  using  only  the  first  term.  Consequently,  we  also  evaluate 

A'  -  MIN  (L). 

®i 

Dj  -  Ej  -  0 

This  procedure  is  equivalent  to  locating  the  radiating  source  which  minimizes  the  square  error  and  then 
comparing  the  minimum  square  error  to  a  threshold. 

To  proceed  with  a  hypothesis  test,  it  is  first  necessary  to  calculate  a  threshold  to  yield  a  desired 
probability  of  false  alarm  equal  to  a.  The  thresholds  T and  f  are  set  so  that 

Pr  {A  >  T\HX)  -  a 

and 

Pr  (A'  >  r\Hx\  -  a. 

In  this  report  we  have  selected  a  »  0.01.  Because  of  this  high  false-alarm  rate,  Monte  Carlo  simula¬ 
tions  can  be  used  to  estimate  T  and  f.  Two  thousand  repetitions  were  run  for  the  case  where  there 
were  n  -  8  array  elements  and  m  -  32  range  cells.  The  signal-to-noise  ratio  (S/N)  was  20  dB  at  the 
output  of  the  antenna  and  the  signal  exhibited  independent  Rayleigh  fluctuations  range  cell  to  range 
cell.  For  each  repetition  the  jammer  location  was  located  randomly  (uniformly  distributed)  within  the 
main  beam  of  the  array.  The  estimated  values  for  the  threshold  were  523  for  f  and  103  for  T. 

PROBABILITY  OF  RESOLUTION -NONCOHERENT  SOURCES 

The  probability  of  resolving  two  Rayleigh  fluctuating  sources  at  a  false-alarm  rate  of  a  -  0.01  was 
found  by  simulation  and  is  shown  in  Figs.  2  and  3  for  the  likelihood  ratio  and  the  one-term  approxima¬ 
tion  of  the  likelihood  ratio  respectively.  Each  figure  shows  the  probability  of  resolution  as  a  function  of 
S/N  for  source  separations  of  0.1,  0.2,  0.3,  0.4,  and  0.5  beamwidths.  After  comparing  Figs.  2  and  3, 
one  concludes  that  the  approximation  to  the  likelihood  ratio  is  only  1  or  2  dB  worse  than  the  likelihood 
ratio  results.  Figure  4  compares  the  likelihood  results  with  those  obtained  [3J  for  an  adaptive  array, 
linear  prediction  least  mean  square  superresolution  algorithm.  The  likelihood  results  correspond  to  a 
0.9'probability  of  resolution.  The  linear  prediction  results  correspond  to  the  appearance  of  a  double 
peak.  Furthermore,  the  number  of  range  cells  is  variable  for  the  linear  prediction  method— enough 
cells  are  used  so  that  the  asymptotic  behavior  is  obtained.  (The  performance  does  not  increase  with  the 
number  of  range  cells  since  the  method  has  an  inherent  noise  bias.  A  nonbiased  eigen-analysis  algo¬ 
rithm  (61  would  have  improved  performance.)  For  sources  separated  by  0.P  beamwidth,  the  likelihood 
approach  requires  approximately  15  dB  less  S/N.  For  larger  separations,  the  difference  is  less. 

PROBABILITY  OF  RESOLUTION-COHERENT  SOURCES 

If  the  two  sources  are  coherent  (for  instance,  one  source  is  just  a  multipath  bounce),  the  signals 
are  not  independent,  but  are  related  by 


PROBABILITY  OF  RESOLUTION 
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S/N  (dB) 

Fig.  2  —  Probability  of  resolving  two  noncoherent  Rayleigh  fluctuating  sources  as  a  function  of 
the  S/N  at  the  output  of  the  array.  Generalized  likelihood  ratio  procedure;  source  separations 
are  0.1,  0.2,  0.3,  0.4,  and  0.3  beamwidth. 


Fig.  3  —  Probability  of  resolving  two  noncoherent  Rayleigh  fluctuating  sources 
as  a  function  of  the  S/N  at  the  output  of  the  array.  Approximation  to  general¬ 
ized  likelihood  ratio  procedure;  source  separations  are  0.1,  0.2,  0.3,  0.4,  and  0.3 
beamwidth. 
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Fig.  4  —  Comparison  of  the  resolution  capability  of  the  generalized 
likelihood  ratio  test  with  a  superresolution  estimating  algorithm 


and 

4>j2  ■  4>]\  +  P  +  (n  -  l)w(sin  9\  -  sin  02)/2.  (8) 

This  complicated  phase  equation  is  required  so  that  the  phase  difference  is  P  at  the  center  of  the  aper¬ 
ture.  For  the  noncoherent  case,  we  minimized  L  [Eq.  (2)]  with  respect  to  4m  +  2  parameters: 
(Aji,  Aj2,  j  “  1.  ....  m),  9\  and  02.  However,  we  were  able  to  solve  for  the  4m  parameters 

in  terms  of  9 x  and  02.  Consequently,  only  a  two-dimensional  search  was  required.  For  the  coherent 

case,  we  must  minimize  L,  Eq.  (2),  with  respect  to  2m  +  4  parameters:  (AJX,  j  -  1 . m), 

p,  p,  0](  and  02-  In  this  case,  one  would  be  able  to  solve  for  the  2m  parameters  in  terms  of  p,  p,  02, 
and  02.  Consequently,  a  four-dimensional  search  would  be  required  and  this  would  require  a  rather 
long  computation  time.  To  avoid  this  long  computation,  we  will  use  the  likelihood  for  noncoherent 
sources  for  the  coherent  source  case.  Obviously,  this  is  a  suboptimal  procedure. 

The  probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  at  a  false  alarm  rate  of 
0.01,  using  the  likelihood  ratio  for  noncoherent  sources,  was  found  by  simulation  and  is  shown  in  Figs. 
5  and  6  for  phase  differences  of  p  -  0°and  p  -  90°.  Figures  7  to  10  show  the  results  for  the  one-term 
approximation  to  the  like  lihood  ratio  for  noncoherent  sources  for  phase  differences  of 
P  —  0°,  p  -  45°,  p  -  90®,  and  p  -  180°.  Comparing  the  results  for  the  two  tests,  one  concludes  that 
the  likelihood  ratio  is  between  2  and  4  dB  better  than  the  one-term  approximation  of  the  likelihood 
ratio. 
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Fig.  5  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  function  of  the 

S/N  at  the  output  of  the  array.  Generalized  likelihood  ratio  procedure;  s  «eoarations  are 

0.2,  0.3,  0.4,  and  0.S  beamwidth.  Phase  difference  is  0°. 


Fig.  6  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  as 
a  function  of  the  S/N  at  the  output  of  the  array.  Generalized  likelihood  ratio 
procedure;  source  separations  are  0.1,  0.2,  0.3,  0.4,  and  0.5  beamwidth.  Phase 
difference  is  90*. 
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S/N  <dB) 

Fig.  7  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  as  a  func¬ 
tion  of  the  S/N  at  the  output  of  the  array.  Approximation  to  generalized  likelihood 
ratio  procedure;  source  separations  are  0.2,  0.3,  0.4,  and  0.S  beamwidth.  Phase 
difference  is  0°. 


Fig.  8  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  as  a  func¬ 
tion  of  the  S/N  at  the  output  of  the  array.  Approximation  to  generalized  likelihood 
ratio  procedure;  source  separations  are  0.1,  0.2,  0.3,  0.4,  and  0.5  beamwidth.  Phase 
difference  is  45*. 
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Fig.  9  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  as  a  func¬ 
tion  of  the  S/N  at  the  output  of  the  array.  Approximation  to  generalized  likelihood 
ratio  procedure;  source  separations  are  0.1,  0.2,  0.3,  0.4,  and  0.3  beamwidth.  Phase 
difference  is  90°. 


10  —  Probability  of  resolving  two  coherent,  Rayleigh  fluctuating  sources  as  a  func¬ 
tion  of  the  S/N  at  the  output  of  the  array.  Approximation  to  generalized  likelihood 
ratio  procedure;  source  separations  are  0.1,  0.2,  0.3,  0.4,  and  0.S  beamwidth.  Phase 
difference  is  180°. 
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SUMMARY 

The  problem  of  resolving  radiating  sources  in  angle  is  formulated  as  a  binary  hypothesis  test.  A 
generalized  likelihood  ratio  test  was  developed  for  noncoherent  sources  and  used  to  resolve  both  non¬ 
coherent  and  coherent  sources.  Two  20-dB  noncoherent  Rayleigh  fluctuating  sources  separated  by  0.1 
beamwidth  can  be  resolved  at  a  resolution  probability  of  0.9  and  at  a  false  alarm  rate  of  0.01  by  using 
samples  from  an  8-element  array  and  32  range  cells.  The  likelihood  method  was  compared  and  found 
to  be  superior  to  a  linear  predictive  least  mean  square  superresolution  estimation  algorithm.  A  one- 
term  approximation  to  the  generalized  likelihood  ratio,  which  is  equivalent  to  the  square  error  residue 
from  fitting  one  source  to  the  data,  is  only  slightly  less  accurate  than  the  likelihood  approach. 
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Appendix 

MAXIMUM  LIKELIHOOD  ESTIMATES 


The  likelihood  function  is  given  by 


L  "  2j  ij  (xu  -  +  CjUt-  DjSt  +  £,7;)2 

+ 1}  ij ( r</  -  cj  w<  -bju>-  ejs<  -  dj  v2- 


If  Hi  is  true,  AJim‘  o  and  hence  Dj  —  E}  —  0.  In  this  case,  Eq.  (Al)  reduces  to 


L  -  Jj  Jj  iXlJ  ~BJW‘  +  C>  U')2 
+  |j  i^av-c^,-^^)2. 


Then  the  maximum  likelihood  estimates  are  found  by  solving  the  equations 

91  91 


BBj  bCj 


Taking  the  derivatives  yields  the  equations 

a  /  n  ft 

-  Jj  -  2 -  Bj w,  +  c>cr,)  +  Jj  -  2t/,(y„  -  c,^  -  *,£/,)  -  o 

and 

|~r  -  Jj  2Ul(XIJ  -BjW,  +  CjU, )  +  Jj —  2Wt(Yff  -  CjW j  -  *,£/,)  -  0. 
Rearranging  terms  and  noting  R'/2  +  i/,2  —  1  yields  the  solution 

Bj  -  Jj  +  DJ  V/», 

Cj  -  jj  W*  -  U,Xu)/n. 

If  Hi  is  true,  the  maximum  likelihood  estimates  are  found  by  solving  the  equations 

AL  AL  AL  AL  0 

BBj  ~  dCj  ~  BDj  “  bEj 
Taking  the  derivatives  yields  the  equations 


jjjj-  Jj  -  2  Wu  ~  BjW,  +  CjU ,  -  DjS,  +  EjT,) 

+  Jj-  lUjiYfj  -CjWj-  BjU,  -  EjS,  -  DjT,)  -  0, 


O.  V.  TRUNK 


|gr- J}2tf,  (*„  -BjWj  +  CjU,  -  DjS,  +  EjT,) 

+  Jj-  2IK,<rtf  -  C,  ^  -  £,S,  -  0,7*,)  -  0. 


+  Jj  -  2T,(Yu  -CjWt-  BjV,  -  EjS ,  -  Z>;r()  -  0, 
jjr-i J  2  w  -  ^  +  C;  £/,  -  DjS,  +  Ej  T,) 

+  Jj  -  2S,(y</  -  Cy  W,  -  fl,t/y  -  £,Sy  -  DjT,)  -  0. 


Making  the  following  substitutions 

F -£(*,$+  V,T,). 
G-T  ~  TO. 

Hj-  ij(^+ 
h  -ifWYu-  U,Xj), 

/,-£<$*„+  r,r„). 
*,  -  2j  ~  W, 

and  rearranging  terms  yields 

nBj  +  DjF  +  EjG  -  Hj, 
nCj  —  DjG  +  EjF  “  !j, 
BjF  --  CjG  +  nDj  *■ 
and 

BjG  +  CjF  +  —  Kj. 


Defining  the  following  vectors  and  matrices 

N  «  ("  0) 


AY_f  ^  <?) 
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allows  £q.  (A2)  to  be  rewritten  as 

NP  +  MQ-  R 
MrP  +  NQ-  Z, 

where  MT  is  the  transpose  of  matrix  M.  Solving  Eq.  (A3)  yields 

P-  (NR-  MZ)/(NN  -  MM7) 

and 

Q  —  (NZ  —  MtR)/(NN  -  MTM). 

Substituting  for  the  matrices  and  vectors  yields 

Bj  -  (nHj  -  FJj  -  GKj)/(n2  -  F2  -  G2), 
Cj  -  (nlj  +  GJj  -  FKj)Kn 1  -  F2  -  G2), 

Dj  -  (nJj  -  FHj  +  GIj)/(n2  -  F2  -  G2), 

and 

Ej  -  ( nKj  -  GHj  -  FIj)/ (n2  -  F2  -  G2). 


13 


